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t^ 1 ■ Abstract 

O ' 

■ There exists an increasing evidence supporting the picture of the Josephson junction (JJ) as a "macroscopic 

quantum system". On the other hand the interpretation of experimental data strongly depends on the assumed 
theoretical model. We analyse the possible states of a Cooper pair box ("charge qubit") for the two types of 
models : two-mode Bose-Hubbard model with its large N aproximations and the many-body description within 
the mean-field approximation (Gross-Pitaevski equation). While the first class of models supports the picture of 
JJ being a quantum subsystem of a single degree of freedom, the second approach yields an essentially classical 
structure of accessible quantum states which, in particular, implies the absence of entanglement for two coupled 
JJ's. The arguments in favor of the mean-field theory are presented and different experimental tests including a 
^-j- ■ new proposal are briefly discussed. 

00 . . 
O ■ 1 Introduction 

O , 

The controversy concerning the classical or quantum character of the Josephson junction (JJ) has a long history. 
In the recent years the picture of JJ as a "macroscopic quantum system" [1] prevails and even more, JJ technology 
is considered to be a very promissing implementation of quantum information processing [2j,||3]. The experimental 
demonstration of "quantum superpositions" of states for a single J J [4J, spectroscopic evidence of coupling between 
two such devices [5] and finally the very recent "measurement of the entanglement of two superconducting qubits 
via state tomography" [6] seem to strongly support this point of view. However, some new theoretical arguments 
in favor of the opposite view - "JJ is essentially a classical system"- appeared also [7]. 

The aim of this paper is to analyse carefully the different mathematical models of JJ using the so-called Cooper 
pair box (CPB) as a concrete example. We begin with the two-mode Bose-Hubbard model and show its equivalence, 
for a large number of Cooper pairs N, to another one - the quantum-phase model with proper boundary conditions. 
These two models essentially support the "naive" picture of quantisation of a nonlinear superconducting circuit 
which leads to a quantum nonlinear oscillator (quantum pendulum model) [8] • Lowest lying two energy eigenstates 
of such a system yield a qubit in the standard approach. 

Then a mean-field many-body theory based on the Gross-Pitaevski equation is applied to a single CPB and two 
coupled CPB's. The obtained structure of the accessible many-particle states is shown to be completely different 
from the predicted by the previous theories. A single CPB cannot support a qubit and the coupling between 
two CPB's introduces classical correlations instead of quantum entanglement. Strictly speaking the amount of 
entanglement can be of the order of 1/N, only. 

To compare both approaches a new Hamiltonian obtained as a large N limit of the two-mode Bose-Hubbard 
model is derived. This Hamiltonian describes a certain nonlinear oscillator model with coherent vectors in a deep 
semiclassical regime directly corresponding to quantum states of the mean-field model. It is argued that only those 
coherents states can be observed in experiments in contrast to unstable eigenvectors of the Hamiltonian. 

Finally, the alternative interpretations of a number of existing experimental tests and the proposal for a new 
one are briefly discussed. 
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2 The theories of Cooper pair box 

We consider a simplified model of JJ assuming that Cooper pairs can be treated as a bosonic gas below the 
critical temperature of Bose-Einstein condensation. We have two electrodes made of a superconducting material 
separated by a thin layer of an insulator which allows for tunneling of Cooper pairs. In order to construct a 
superconducting qubit one tries to suppress the tunelling of many Cooper pairs by using the circuit consisting of a 
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small superconducting island "1" connected via J J to a large superconducting reservoir "2". Coulomb interaction 
between Cooper pairs in a small electrode become important and must be taken into account in the Hamiltonian. 

2.1 Two-mode Bose-Hubbard model 

The annihilation and creation operators Si, a} and a 2 , a\ correspond to the ground states of a boson (Cooper pair) 
in separated electrodes "1" and "2", respectively. Coulomb interaction between Cooper pairs in a small electrode 
is modelled by the quadratic term in the Hamiltonian below. The second term contains the potentials of the both 
electrodes and the term proportional to K describes the tunneling of Cooper pairs 

H 2 ,nod = EciaXai) 2 + ~(JJx<i\a\ + J7 2 a|a 2 ) - K{a x a\ + a\a 2 ) . (1) 

The total number of bosons if fixed and equal to N, hence we should consider the restriction H' 2mod of the Hamil- 
tonian |T]) to the N + 1-dimensional Hilbert subspace of physical states spanned by the orthonormal basis 

|fc) = WN-k)\]-^(a\) k (ai) {N - k) \vac) . (2) 

H' 2mod is defined as follows (we omit an irrelevant constant and use the convention |JV + 1) = | — 1) = 0) 

KnoAk) = [Ec(k - n x ) 2 ] \k) - K^k{N-k){\k + 1) + |fc - 1)) , k = 0, 1, N. (3) 

where U = U 2 — U\ and h\ = U/AEc is an averaged number of Cooper pairs on the island corresponding to a 
stationary reference state. As we are interested in the low lying energy eigenstates we can replace in the second 
term K^Jk{N — k) by K^ni{N — hi) = Ej, Therefore, the final form, which can be called the Bose-Hubbard 
Hamiltonian, reads 

H B - H \k)= [E c (k-hi) 2 ]\k)-Ej(\k + l) + \k-l)) , fc = 0,l,...,iV. (4) 

with the convention \N + 1) = |0), | — 1) = \N). This Hamiltonian can be interpreted as describing a particle with a 
discretized position A: in a harmonic potential located around k = hi and a hopping term modelling its free motion. 

2.2 From two-mode to quantum phase model 

The standard picture of the CPB can be obtained from the Hamiltonian ([I]) treating the operators hi , 02 as classical 
variables ai , a 2 and using the following parametrization 

ni = 0*01, n 2 = a 2 a 2 , ai = y/nle 101 , a 2 — ^/n^e 1 ® 2 , 6 2 — 81 = 9, n\ + n 2 = N, (5) 

which, under the condition |m — n%\ << fix, yields the Hamiltonian (up to an irrelevant constant) 

H = E c (ni - fn) 2 - Ej cos6 . (6) 

The next step in this rather "naive" approach is to "requantize" the variables m, assuming the following canonical 
commutation relations 

[h 1 ,6]=i, (7) 
and to consider the properties of the "quantum nonlinear oscillator" with the Hamiltonian 

H = E C £ 2 -Ej cos§ , (8) 

where ^ = hi — fi\ is a "momentum" operator conjugate to 0. The two lowest lying eigenstates of H, which are 
automatically orthogonal, are supposed to provide the superconducting implementation of a qubit. 

A more subtle and mathematically sound approach starts with the Hamiltonian ([4]) and treats the basis \k) as a 
subset \k) = (l/v / 27r)e I (' £- " 1 ) 6 '; k = 0, 1, ...N of the (rotated by e~ mid ) Fourier basis in the quantum rotator Hilbert 
space L 2 [— it, n]. Assume that hi is a natural number. Then on the subspace of wave functions f(0) spanned by 
these particular basis elements 

jV N-ni 

f{0) = Y, Ckei{k ~ ni)(> = E c fc+fil e lfce ,c fc GC (9) 

k—0 k——ni 
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the action of the operator Hb-h is equivalent to the action of the following differential operator 



H qph f(e) = -E c ^f(0) -Ej cos 9f(6) (10) 

with the periodic boundary condition f(ii) = /(— 7r). Then we treat the operator U qp h extended to the whole Hilbert 
space L 2 [— tt, tt] as a good approximation to the original Hamiltonian |(4|). This is, indeed, a reasonable procedure 
because N >> hi and hi is at least of the order of 10 s . Hence, the restriction to a finite sum in J9]) becomes 
irrelevant at least for the interesting low energy regime. We can continuously change the potential difference U to 
produce a continuous change of the parameter hi = U/4Ec beyond the natural numbers. This disagrees with the 
necessary periodicity of the functions J9]) from the domain of the Hamiltonian H qp h . Such a generic situation leads 
to the slightly modified Hamiltonian (a = hi — integer part[n\]) 
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H qph f(6) = -E c - ia\ f{9) - Ej cos Of {6) . (11) 
2.3 Mean-field many particle model 

The fundamental property of the Bose-Einstein condensate is the form of its time-dependent N-boson wave function 
ty{t) (i.e. the true quantum state of the system) which is approximately given by a tensor product of the identical 
single-boson wave functions ip(t) |9J 

*(*) = ®V»(*)- (12) 

N 

One should notice that the set of condensate's wave functions does not form a linear subspace in the Hilbert space 
of the N-particle system. 

There are numerous arguments, including theorems in mathematical physics [10J,[11J that for weakly interacting 
Bose gas the initial product state evolves approximately into the product state (fT2"ll such that the single-particle wave 
function tp(t) satisfies the Gross-Pitaevski equation which is a nonlinear Schrodinger equation with a cubic term 
modelling boson-boson interaction |13| . It implies immediately that the wave function ip(t) does not correspond 
to a quantum state of any subsystem and the superposition principle for ip{t) is valid only approximatively, due 
to nonlinearity. Hence, ^p(t) is rather a classical object: the order parameter labelling different thermodynamical 
phases of the system below the critical temperature. 

For our model the spatial dependence of the wave function ip is irrelevant and therefore we represent ip by the two 
complex amplitudes -0 = (ipi,^) corresponding to the electrodes "1" and "2", respectively. They satisfy a discrete 
version of the Gross-Pitaevski equation which, for the convenience, is formulated in terms of the renormalized wave 
function </> = y/Nip [12] 

0i= -i(Ui + g\cj)i\ 2 )(j)i - iK4> 2 

02 = -iU 2 h - iKfa (13) 

with 1 4>\ | 2 + 1 4>2 1 2 = N. The term proportional to g describes Coulomb interaction of Cooper pairs in the first 
electrode in the mean-field approximation. 
Introducing the following parametrization 

<f>j = Qj + *Pi = V^^ 3 . 3 = !> 2 ( 14 ) 

one can easily check that the equations lfl3|) can be written as the Hamiltonian equations for the system of 2-degrees 
of freedom governed by the classical Hamiltonian 

h 2 (qi,q2,Pl,P2) = — 2 U2 (Pi + g?) + ~ (p\ + qlf - K( Pl p 2 + gl g 2 ) + ^(p 2 + q\ + p 2 2 + q\) . (15) 

The system of above is integrable because of the existence of the constant of motion ||0|| 2 = {p\ + q\ + p\ + q 2 ) = N. 
We can restrict ourselves to a decoupled, single degree of freedom subsystem parametrized by 9 — 2 — Q\ and n\ 
(fT4|) with the Hamiltonian function 

9, U, 



hi(9,m) = 2.(m ) 2 - KJriJN ~ ni )cos9 (16) 

4 g 
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where U = U2 — U±. Again, in the regime \n\ — n\\ << -hi the Hamiltonian (fT6| coincides with the semiclassical 
limit J6]) of the 2-mode Hamiltonian if g/4 — Ec, n\ = U/g and Ky/fi\(N — WTj = Ej, We can solve the 
corresponding Hamiltonian equation to obtain a trajectory (0(f), £(f)) with the "momentum" £(f) = n\[t) — fi\. 
Taking into account that £(t) is at most of the order y/Ej/ Ec and eliminating the irrelevant common phase factor 
we can use the following approximation for <f> = (<pi, $2) 



i8(t) 



ni + 



- ni)m 



Mt)= y/N-fi!- W) ~ -£(*)( 2^ 



y/N-fi! 



(17) 



Now, the normalized to one single boson wave function ip — N 1 / 2 <f> can be written as the superposition of two 
orthonormal vectors rotated by the unitary W 

m = W(9(t)) [ Xo + £(t) (2y/(N - ni)fn) 

where 

W(9) =diag(e- 4e ,l) , Xo-^ 1/2 (V^T, y/N-n^, Xl = N-^^y/W 



711,- 



(18) 



(19) 



One should notice that the amplitude of the vector xi is scaled down by the factor (2y/ (N — ni)ni) _1 << 1 and 
hence the normalization of ip(t) is preserved up to the higher order terms. 

We can reconstruct now the quantum state \I/(t) describing the condensate in the mean-field approximation 
keeping the lowest order correction of the order 1 / yfWy 



*(t) - W(0(t)) [$ + C(*)Viv(2V(iV - ni)^ 



$1 



where 



W(0) = (g) W{0) , $ = (g) Xo , $1 = 

JV iV 



1 



A? 

E 



(20) 



(21) 



® Xo ® Xi ®Xo ® 
• 7=1 j-th 

One should notice that the state $1 does not describe the condensate but the states of the form $0 + e$i do, 
up to the order ~ e 2 . It is rather clear that the obtained 2-dimensional manifold of the available states does not 
contain the manifold spanned by the qubit states of the form a\0) + f3\l) with \a\ 2 + |/3| 2 = 1. Another equivalent 
interpretation of the product states (fT2l) in terms of coherent vectors will be presented in Section 4. 



3 Phenomenology of the mean-field model 

The experiments performed on the different types of "JJ qubits" are usually analysed in terms of the standard 
quantum-phase model. The question arises whether the phenomenology of JJ can be also consistently explained 
within the many-particle mean-field model. 



3.1 Single Cooper pair box 

The quantum state of the controlled single CPB is described by the formula f20|)(|2T1) with 0(f), £(t) being the 
solutions of the classical Hamiltonian equations obtained for the Hamiltonian (fl6|) with a generally time-dependent 
U(t) given by 

U(t) = U + u(t) . (22) 

Here U is a constant potential including the constant external voltage and the interaction of Cooper pairs with 
the background positive charges. The term u(t) describes external control usually in the form of pulses or periodic 
perturbations. Therefore the effective classical Hamiltonian can be written as 

h{9, £; f) = Ec? - tt(t)f - Ej cos 6 . (23) 

Starting with the stationary state ^(0, 0) and applying a relatively weak pulse u(t) we can excite the small oscillations 
of 9(t) and £(f) with the plasma frequency ujq ~ \J2EqEj . From the form of the quantum state ifTTj) it follows 
that the excess charge on the electrode " 1 " is equal to 2e£ (f ) and hence oscillates with the same frequency emitting 
an observable radiation. A stronger pulse increases the " kinetic energy" E c ? above the "potential barier" 2Ej 
leading to a different regime of the dynamics. Applying a periodic field one can create a variety of resonanse 
phenomena and even chaotic behavior. 
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3.2 Two coupled CPB's 



The coupling of two CPB's can be realized by the electromagnetic interaction between the oscillating charges of the 
electrode "1" of the first CPB and the electrode "1"' of the second CPB which in the Gross-Pitaevski equations 
of below is described by the terms proportional to G. The coupled equations for the two wave functions of the 
condensate cp, </>' read 

k= -i{Ui+g\<i>i\ 2 ^CW^^-iK^ 
4>2 = -iU 2 4>2 - iKcf>i 

<f>'i = -iiUi+g'W^ + G^M-iK'fa 

4>' 2 = -iU&'z - iK'4{ . (24) 

By the procedure, completely analogical to a single CPB case, one can obtain the effective Hamiltonian for two 
degrees of freedom described by two angles 8, 9' and momenta £, £' 

h(9, 9', £, O = E C C 2 - Ej cos 9 + E' c i' 2 - E'j cos 9' + . (25) 

Here, analogically to a single CPB case: U = U2 — U1 , U' = U 2 — U[, £ = n\ — n\, £' = n'i — n[, Ec = <?/4, E' c = g'/A, 
Ej = K^(N — ni)ni, E'j = K'y/ (N' — n'-^n^ and n\,fi\ are solutions of the equations U = gfi\ + Gn[,U' = 
g'fi! x + Gn\ . 

Solving the Hamiltonian equations for the coupled nonlinear oscillators with possible external control we obtain a 
classical trajectory {9(t),9'(t), £(£),£'(*)} which defines the evolution of a quantum state of the coupled two CPB's. 
This state is a tensor product of quantum states for separate CPB's, each of them given by ([20|) 

€,O = *(^0®*(^O • (26) 

Therefore, it is not possible to produce entangled states within this model. A small amount of entanglement of the 
order 1/N, which is allowed according to [7], is well-placed within the accuracy limits of the mean-field approach. 
Obviously, the quantum states of the two JJ's are classically correlated by the coupling of the classical parameters 
6,£ and 6',?. 



4 Comparison of the models 

We have obtained two physically different pictures of the CPB. According to the first one we can observe and 
manipulate lowest lying states and their superpositions for the well-defined quantum system equivalent to the 
quantum pendulum. The second one implies that the accessible states form a nonlinear manifold parametrized by 
the essentially classical variables satisfying classical equations of motion. To compare both approaches we represent 
the product state (fT2|) in terms of the basis introduced in ([2]) 



TV 1 



1 n i -t 
1 al, 

N 2 



N 



N 

E 

fc=0 



\NJ V TV 



N-k 



\k) 



(27) 



Replacing now the square root of the binomial distribution in l|27|) by the square root of the Poisson one (what is 
consistent with the assumption N >> n\) we can treat the states \k) as the eigenstates of wb for a certain fictitious 
nonlinear oscillator and the product state ([271 becomes the coherent vector 



I*) 



E 

fc=0 



-i0\ 



where 



Using now the identities 



(28) 



(29) 



a\d 2 \k) = VN -kP\k) 



a\ai\k) 



lb\k) 



(30) 

and the fact that the relevant states are spanned by the vectors \k) with k ~ hi we can approximate the Hamiltonian 
(U) by a new one 

H 2mod ~H b = E c {Pb - hi) 2 - -§-(S + &t) . (31) 
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Here again we obtain the Hamiltonian which in the semiclassical limit yields the same equations of motion as 
J6]) . In contrast to the previous approximative Hamiltonians (J4j) ©) , where the large values of the quantum number 
k are hidden by applying the shift k i— » k — n\, here it is clear that the lowest eigenstates of Hb are placed in a 
deep semiclassical regime corresponding to the quantum numbers (eigenvalues of b'b) of the order of fi\ ~ 10 8 . 
Therefore, we should not expect to observe in experiments the true eigenstates of the Hamiltonian H& but rather 
localized "classical" states described by the coherent vectors |28|) evolving (approximately) according to the classical 
Hamiltonian equations (see the next Section). The striking difference between the Hamiltonian eigenstates and the 
coherent vectors concerns the fluctuations of the Cooper pair number n\ and the phase difference 9. For the 
eigenstates both fluctuations are of the order 0(1) while for the coherent states we have normal fluctuations of n\ 
of the order ^/n{ and for the phase of the order of . The later behavior is consistent with the fact that 6 is 
an essentially classical variable - the order parameter for the Bose-Einstein condensation of Cooper pairs. 

5 Stability of classical states 

The reasons for a "spontaneous transition to the classical world" in the regime of large quantum numbers have 
been extensively discussed since the famous "Schrodinger Cat Paradox" [14]. Namely, the quantum states of large 
systems can be observed only if they are sufficiently stable with respect to the interaction with an environment. For 
large enough systems the stable states are described rather by the strongly localized symmetry-breaking solutions of 
the nonlinear mean-field type equations (e.g. Hartree/Hartree-Fock eqs., Born-Oppenheimer approximation) than 
by the eigenvectors of the exact Hamiltonians. This is for instance the case of fixed molecular structure for large 
enough molecules (e.g. famous chiral molecules) or deformed nuclei. For the similar reasons nonlinear mean-field 
Gross-Pitaevski equation should give a proper descriptions of the observable states of N Cooper pairs in contrast 
to the eigenvectors of the isolated system Hamiltonians |T]) , ([3]) , |4j , © , or (|31|) . 

The complete analysis of the stability of the models given by (U) and (f3Tj) will be presented in the forthcomming 
paper. Here we discuss briefly the model lf3lj) assuming that the main source of instability is the transition of bosons 
from the condensate to other "normal" modes and back. This effect can be modelled by adding the dissipative part 
to the equation of motion for the density matrix of the oscillator 

j t P = -i[H b ,p]+L(p) , (32) 

L(P) = b f ] + [b, 0]) + \s{pp, b] + p,pb]) . (33) 

A simple measure of stability of a pure state \<j>) is given by the initial decay rate of fidelity due to the dissipative 
part of the dynamics 

r(|$)=-<0|L(|0W|)|0 =(7 + ^X(#%>H(#l0}| 2 )+^ (34) 

The decay rate is minimal and equal to 6 for the states satisfying the eigenvector condition b\<p) = a\4>) i.e. for 
the coherent vectors f29|) . The Fock vectors decay with the rates r(|fc)) = (7 + S)k + S what means that the 
life-time of the low lying states of (|3T|) which are supposed to support "superconducting qubit" is about n\ ~ 10 8 
times shorter than the life-time of classical coherent states. Therefore, the states observed in experiments should 
be described by mixtures of coherent vectors. 

6 Comparison with experiments 

We discuss now briefly several, already performed experimental tests aimed to justify the picture of superconducting 
qubits. 

Coherent oscillations. The charge oscillations for CPB have been confirmed for the first time in [4]. The mean- 
field model predicts small oscillations with the plasma frequency loq = \/2EcEj while for the charge qubit model 
the frequency is a more complicated function of Ec,Ej and the applied voltage U [3]. It is not clear whether this 
difference can be tested experimentally as it seems that no independent measurements of Ec and Ej are available 
and moreover JJ's parameters are subject to time-dependent fluctuations [2]. 

Transmissivity of CPB. The interesting experiment [15] shows the existence of two distinct states of CPB 
correlated with two values of transmissivity for CPB connected to a single electron transistor. However, classical 
oscillator shows also a "two state structure" described by the (averaged in time) probability distribution for its 
canonical variables - p(x) = ~(1 — a; 2 ) -1 / 2 ( amplitude = 1) - with two infinite maxima at ±1. 
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Particle number fluctuations. It seems that the large fluctuations (<~ \/ni) of the Cooper pair number predicted 
by the mean-field theory contradict the experimentally proved precise control of n\ [16]. However, this experiment 
has been performed with the superconducting-normal junction. Therefore, Cooper pairs exist on the superconduct- 
ing island only and the normal electrode acts as a heat baths. It follows, that in this experimental setting the state 
of Cooper pairs is given by the Gibbs density matrix 

p=Z- 1 exp{-^(a t 1 a 1 -n 1 ) 2 } (35) 

which at low temperatures (kT < Ec) exhibits, indeed, very small fluctuations of a\a\. 

Rabi oscillations, Ramsey fringes. The observation of Rabi oscillations and the associated Ramsey fringes are 
easily and elegantly described by JJ-qubit models \S\. On the other hand , as shown in [17j,[18j,[19|, the experimental 
data can be attributed to classical nonlinear dynamics and statistics as well. 

Entanglement via state tomography. The experimental data for two coupled current-biased JJ's has been analysed 
in terms of the "state tomography of two phase qubits" [6]. Here again the analysis is strongly model-dependent 
and does not exclude the alternative interpretation in terms of coupled nonlinear oscillators. 



7 A new test of quantumness 

The presented above experimental tests, although quite convincing, are not ultimate because of their strong model- 
dependence. Model independent Bell inequalities are often proposed as an ultimate test of quantumness. In the 
case of coupled JJ's the corresponding experiments are still to be done. Unfortunately, even a positive result - 
violation of Bell inequalities - can be challenged because of the strong presence of the "locality loophole" in this 
setting. Namely, in contrast to experiments involving photons we are not able to separate JJ's after the act of 
interaction. Another interesting possibility is to test Bell inequalities for multi-time correlations in a single system 
as proposed in |20J. 

Recently, a new model- independent test of quantumness applicable to a single system has been proposed [21] . 
Assume that the CPB is indeed a two-level system. Then we can choose the following observables A and B given 
in a fixed basis by the following matrices: 

/0.724 0.249 \ „ A \ 

^0.249 0.0854,1 ' \0 0.309,1 ( } 

It is easy to check that for an arbitrary state ip the following inequalities hold 

0<(il>\A\il>)<(1>\B\il>). (37) 

On the other hand for a specially choosen state 4> — 920) W6 °^ tam ' < ( t'\^ 2 \ ( t > > ~ < 4>\A 2 \(t) >= —0.0590. This 

clearly demonstrates the quantum nature of the qubit because, for any classical model representing the observables 
A, B by the functions on a certain "phase space" and states by the probability measures, if the inequality < (A) < 
(B) holds for all probability measures then always < (A 2 ) < (B 2 ). 



8 Conclusions 

The large body of experimental results which can be easily and elegantly interpreted in terms of superconducting 
qubits suggests to accept the quantumness of Josephson junctions. On the other hand, as discussed above, the idea 
of "macroscopic quantum system" applied to N-particle systems with N of the order of 10 8 deeply contradicts our 
understanding of their quantum and statistical properties. In particular, it is generally believed and experimentally 
confirmed in chemistry, nuclear physics , quantum optics and physics of critical phenomena that already for much 
lower particle numbers or quantum numbers, the eigenvectors of the Hamiltonians are highly unstable with respect 
to the interaction with an enviroment. As a consequence, the observed quantum states are perfectly described by 
the stable localised symmetry breaking solutions of the non-linear mean-field type equations. These solutions satisfy 
classical equations of motion with additional damping terms. The system of 2N electrons in a superconducting state 
seems to be ideally suited for the mean-field description in terms of the Gross-Pitaevski or rather the Ginzburg- 
Landau equations. Therefore, it is still necessary to discuss alternative interpretations of existing experiments and 
search for new ones until clear evidence to quantum or classical character of JJ's becomes available. The proposal 
of such an ultimate experiment is discussed in this paper. 
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